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Abstract 

It is proposed to create materials with a desired refraction coefficient in 
a bounded domain D C R'^ by embedding many small balls with constant 
refraction coefficients into a given material. The number of small balls 
per unit volume around every point x £ D, i.e., their density distribution, 
is calculated, as well as the constant refraction coefficients in these balls. 
Embedding into D small balls with these refraction coefficients according 
to the calculated density distribution creates in D a material with a desired 
refraction coefficient. 
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1 Introduction 

In [I]- [3] it was proposed to create material with a desired refraction coefficient 
by embedding into a given material small particles with suitably chosen bound- 
ary impedances. It was proved that any desired refraction coefficient n^(x), 
£jn^(x) > 0, can be created in such a way in an arbitrary given bounded domain 
D C M^. Preparing small particle with a prescribed large boundary impedance 
may be a technologically challenging problem. By this reason we propose in this 
paper a new method for creating materials with a desired refraction coefficient. 
We hope that this method may be easier to implement practically. 

This method for creating materials with a desired refraction coefficient n'^{x) 
consists of embedding into a given material small particles (balls) with a suitably 
chosen density distribution of the embedded particles and a suitably chosen 
constant refraction coefficients of each of the embedded particles. No boundary 
impedances are necessary to use in this method. Therefore, one hopes that the 
new method may be easier to implement in practice. 
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The density of the distribution of the embedded particles and their constant 
refraction coefficients are calculated given the desired refraction coefficient n'^(x) 
and the refraction coefficient nQ{x) of the original material in D. 

Let us formulate the problem precisely. Assume that a bounded domain D C 
M.^ is filled with a material with known refraction coefficient rig (a;), '^tiq (x) > 0, 
ng(x) — 1 in D' :— M.^ \ D, nQ{x) is piecewise-continuous. Throughout this 
paper by piecewise-continuous function we mean a bounded function with the 
set of discontinuities of Lebesgue measure zero in R'^, and do not repeat this. 

The waves satisfy the equation: 

Lquo ■■= [V^ + k'^nl{x)]uQ = m R^, k = const > 0, (1) 

uo = e*'^-"-^ + V. (2) 
Here v is the scattered field, satisfying the radiation condition: 

Vr — ikv — o{r^^) r := \x\ oo, (3) 

where a € is the direction of the incident plane wave, and S"^ is the unit 
sphere in M.^. It is proved in [l] that this scattering problem under the stated 
assumptions on nl{x) has a unique solution, and the function G{x, y), satisfying 
the equation 

LoGix,y) = -6ix-y) m M^, (4) 

and the radiation condition ([3]), does exist and is unique. 
One can write 

io = + fc^ - qo{x), 

where 

qo{x) := qo{x, k) := k^ - k'^nl{x), qq^x) =0 x £ D' . (5) 

Let be a desired refraction coefficient in D. We assume that is 

piecewise-continuous, 3ri^(a:) > 0, and n^(x) = 1, a; G D' . 

We wish to create material with the refraction coefficient in D by 

embedding into D many small non-intersecting balls f < m < M, of radius 
a, centered at the points Xm G with constant refraction coefficients in 

Smallness of the particles means that ka « 1. 

Let A C D be any subdomain of D. We assume that the number of small 
particles, embedded in A, is given by the formula: 

A/-(A) ^ 1 = KT' / N{x)dx[l + o(f )] a ^ 0, (6) 

where N{x) > is a piecewise-continuous function in D, and 

K := W/3 
is the volume of a ball of radius a. 
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Formula ^ gives the density distribution of the centers of the embedded 
small balls in D. The total number of these balls tends to infinity as 0{V~^) = 
0{a~^) when a — > 0. 

We assume that the total volume V{D) of the embedded particles (balls) is 
not greater than \D\, where |D| is the volume of D, i.e., 



V{D) = VaMiD)= [ N{x)dx[l+oil)]<\D\, 

JD 



This means physically that although N(x) can be large in some subdomains of 
D, its average over D is not greater than 1. 

The scattering problem in the case of the embedded into D particles is: 



M 



(Lo + ^ n2^Xm)t/ = m R\ (7) 



m — 1 



where Xm is the characteristic function of the ball Bm, i.e., Xm = 1 in Bm-, 
Xm = in R3 \ Bm, and 

U = uo + V, (8) 

where V satisfies the radiation condition uq solves the scattering problem 
in the absence of the embedded particles, i.e., when M — 0, and 

Here v^{x) is some piecewise-continuous function in D, Qv'^^x) > 0. 

The solution U{x) — Ua{x) to problem Q"® depends on the parameter a, 
and the number M of the embedded particles depends also on a, 

M^O{V-^)^0{a-^), 

so M ^ oo at a prescribed rate as a 0. We are interested in the limiting 
behavior of U{x) — Ua{x) as a ^ 0. Our basic result, Theorem 1, below, says 
that the limit 

lim Ua{x) := Ue{x), (9) 

does exist and satisfies the integral equation p7p in Theorem 1. 
From Q-dHl) one gets 

U{x) = uoix) + fc2 V / G{x, y)U{y)dy. (10) 

This integral equation we rewrite as 
j\/ » 

Uix) =uo(x)+ey2nl^ Gix,y)dyU{x„^)[l + o{l)] a ^ 0. (11) 
T Jb^ 
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Here the continuity of U in i3,„, 1 < m < M, was used. This continuity implies 

Uiy)^Uix„,)[l + o{l)] a^O; y e B^. 

The function U is twice differentiable in R^, as follows from pT|) . so it is con- 
tinuous in D. 

Wc need three lemmas. 

Lemma 1. The following relations hold: 

lim \x -y\G{x,y) ^ (12) 
sup |x-y||G(x,2/)l < c. (13) 

[x-j/|>0 

By c > we denote various estimation constants. 
Proof of Lemma 1 is given in Section 2. 
Lemma 2. The following relations hold: 

[ I. = V.I. -.„!-'. i.-.J><.. 

\x - yl^'^dy = 2Tr[a'^ - — — |^^), \x - x^l < a. 

\y—Xm\<a ^ 



(14) 



Proof of Lemma 2 consists of a direct routine calculation and is therefore 
omitted. The result of Lemma 2 is known from the potential theory. 

Lemma 3. // / is piecewise- continuous in D and the points Xm are dis- 
tributed in D by formula then the following limit exists: 

M „ 
lim Va V f{x„,) = / f{x)N{x)dx, (15) 

where N{x) is defined in ([S]). 
This Lemma is proved in [T]. 
One has: 

G{x,y)dy^VaGix,x„,)[l + o{l)l a -> 0. (16) 

From (fTTI) . (fTBjl and our basic result follows: 
Theorem 1. There exists the limit ([9]) and 

Ue{x) = uo{x) + I G{x, y)N{yy{y)ue[y)dy. (17) 

JD 

Physically the limiting field Ue is interpreted as the effective (self-consistent) 
field in D. 
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Corollary 1. The functions U{x) and Ue{x) are twice differentiablc in R'^. 
The function Ue(x) solves the equation: 

Lueix)=0, L:^ Lo + k^N{x)i'^{x), (18) 

so 

ri^ix) =nl{x) + N{x)v^{x). (19) 

To prove this Corollary one applies the operator Lq to equation (fT7|) and 
uses equation (|4]). 

Conclusion: To construct a material with a desired refraction coefficient 
n^(x) one embeds small balls with radius a, ceneterd at the points Xm, 1 < m < 
M, distributed by formula and chooses N(x) andv^ix) so that relation (jl9p 
holds. 

The choice of N{x) and v'^{x) is therefore non-unique, because the relation 
(|19p can be satisfied by infinitely many ways. For example, one may fix N{x) > 
in D and then choose 

2( . ^ n'^ix)-nl{x) 
^ ' N{x) 

If n'^{x) = in a subdomain A C -D, then one can take N(x) —0 in A. 
In Section 2 proof of Lemma 1 is given. 

2 Proofs 

Proof of Lemma 1. We start with the equation: 

G{x,y) = g{x,y) - I g{x, z)qo{z)G{z,y)dz := g - TG, (20) 
Jd 

where go is defined in (O, and 

pik\x—y\ 

g{x,y)^— (21) 

4:Tr\x — y\ 

Equation ([^0]) is of Fredholm-typc in the space X of functions ip{x,y) of the 
form ^{x,y) = where (j){x,y) is a continuous function of its arguments, 

and the norm in X is defined as ||V'|| = sup^. j,gR3 (|x — y\\ip{x,y)\). 

We have ||g|| = j^. The homogeneous equation ([^0)1 has only the trivial 
solution (see [1]), so the operator (J + T)^^ is bounded in X. Therefore, ||G|| < 
c\\g\\ = This implies estimate ([T3l) . 

To prove ([E]), let us multiply ^ by \x - y\ and let \x - y\ -> 0. One 
has lim|j;_j,|^o .9 — 37- The integral TG is bounded for all x,y G D, so 
lim|^_y|_o(|a; - y\TG) = 0. Thus, relation IH]) follows. 

Lemma 1 is proved. □ 
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